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MEIXNER MATRIX ENSEMBLES 

WLODZIMIERZ BRYC AND GERARD LETAC 



Abstract. We construct a family of matrix ensembles that fits Anshelevich's regression 
postulates for "Meixner laws on matrices" , namely the distribution with an invariant 
property of X when E(X 2 |X + Y) = o(X + Y) 2 + 6(X + Y) + cl„ where X and Y are 
iid on symmetric matrices of order n. We show that the Laplace transform of a general 
n x n Meixner matrix ensemble satisfies a system of partial differential equations which 
is explicitly solvable for n = 2. We rely on these solutions to identify the six types of 
2x2 Meixner matrix ensembles. 



^ ■ 1. Introduction 

The classical Meixner laws are the one dimensional binomial, Poisson, negative bino- 
mial, gamma, normal, and hyperbolic Meixner laws that include hyperbolic secant law. 
These laws make their appearance as the orthogonality measures of a certain family of 
orthogonal polynomials [32J, as the laws characterized by quadratic regression property 
[25], and as the exponential families corresponding to the quadratic variance functions 
[33], see also [19]. The Laha-Lukacs [23l [25] characterization is the description of the 
possible distributions of a square-integrable random variable X such that there exist real 
numbers a, 6, c satisfying E>(X 2 \S) = aS 2 + bS + c where S = X + Y, where Y ~ X and 
where Y is independent of X. Another convenient form of this condition is to postulate 
the existence of real numbers A, B,C such that 

(1.1) E((X-Y) 2 \S) = AS 2 + BS + C 

(use (X - Y) 2 = 2X 2 + 2Y 2 - S 2 to see that (A, B, C) = (4a - 1, 46, 4c)). 

In 2006 M. Anshelevich [U pages 22 and 25] observes that (11. ip makes sense when the 
real random variable X is replaced by a random square matrix X. He then raises the 
question of defining analogues of Meixner distributions on matrices. In other words, he 
asks for the matrix version of Laha-Lukacs [25] result. This question on random matrices 
is so general that it is suitable to restrict this problem to random variables X valued in 
the space H n l of all symmetric matrices and such that such that UXU* ~ X for matrices 
U in the orthogonal group O(n) = K, n ^. On the other hand, it is natural to extend this 
simpler framework to the space H n 2 of Hermitian complex matrices for random variables 
X such that l7X[/* ~ X for matrices U in the unitary group U(n) = lC n ^ and even to 
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the space H n 4 of Hermitian-quaternionic matrices, again for random variables X such 
that U~X.II* ~ X for matrices U in the symplectic group Sp(n) = K. nj ±. There are also 
connections with free probability (see [9] and Section [6731 below), and connections with 
Jordan algebras (see Section [6741 below) . 

A probability law /i on H nii g (where = 1, 2 or 4), or a random variable X with law 
//, is said to be rotation invariant if the law of UXU* is also \x for all U in K n ^- In such 
case, following a physicists tradition we shall call [i an ensemble. 

Definition 1.1. We will say that the probability law \x on H nji g is a Meixner ensemble 
with parameters A, B, C G R if /J is rotation invariant, the second moments exist, and \x 
has the following property: if X, Y are independent with the same law /x and S = X + Y, 
then 

(1.2) E((X- Y) 2 |S) = AS 2 + BS + CI n . 

This is a different concept of a Meixner matrix ensemble than the one introduced in 
[IB*] under the name "multivariate Meixner classes of invariant distributions of random 
matrices". In particular, we note that the oldest known ensemble, namely the Wishart 
distribution with mean proportional to I n is not a Meixner ensemble in the sense of 
Definition ll.il since if X, Y are i.i.d. with Wishart distribution, then there exist two real 
numbers A, B such that 

(1.3) E((X-Y) 2 |S) = AS 2 + 5StrS, 

which is quadratic in S but not of the desired form (ll.2p . For a proof of (jl.3p . see [271 
page 582]. 

To complete this introduction, let us make a number of simple remarks about the 
Meixner ensembles. It is easy to check that if the law of X is a Meixner ensemble with 
parameters A,B,C, then an affine transformation X = aX + t with real a,t is also 
Meixner, with parameters 

A = A, B = aB — 4At, C = a 2 C + 4At 2 - 2Bat. 

In particular, since passing to —X changes only the sign of B, without loss of generality 
we will consider only B > 0. Suppose now that E(X) and E(X 2 ) are given; in this case 
the constants A,B,C satisfy an equation with coefficients that depend on E(X), E(X 2 ). 
The equation is obtained from taking the expected values of both sides of (11.20 . Finally, 
rotation invariance implies that both E(X) and E(X 2 ) are constant multiples of I n , so in 
the non-degenerate case E(X 2 ) ^ we can normalize the matrices so that E(X) = and 
E(X 2 ) = l n , and express relation (ll.2p in terms of two real constants a, b as 

(1.4) E((X - Y) 2 |S) = — ^- (l n + bS + aS 2 ) . 

In this normalized setting one can classify the Meixner ensembles into 

(i) the elliptic case b 2 > 4a, which in Laha-Lukacs theorem on R corresponds to 
negative binomial (a > 0), Poisson (a = 0), or binomial (a < 0) laws. 

(ii) the parabolic case b 2 = 4a, which corresponds to gamma law (a > 0) or Gaussian 
law (a = 0) 
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(iii) the hyperbolic case b 2 < 4a, which corresponds to the family of hyperbolic 
Meixner laws that include hyperbolic secant law. 

Of course, for non-degenerate ensembles (11.21) and (11.41) are equivalent. If (11. 2p holds 
with constants A ^ 1,B,C, then the standardization X of X with mean fj,I n and variance 
Var(X) := E(X 2 ) - (E(X)) 2 = a 2 I n satisfies (OD with 

/-i r\ A B + AAfi 

(1.5) c 



2(1 - A) ' 2a(l -A)' 

What we call trivial Meixner ensembles are random multiples of the identity, X = £I n 
where £ is R- valued with one of the six Meixner laws. The aim of the paper is the study 
of non-trivial Meixner ensembles for all n. We will be able to find all of them for n = 2, 
we will give examples of Meixner ensembles for n > 3 and we will write a system of n 
linear partial differential equations (PDEs) of the second order satisfied by the Laplace 
transform of a Meixner ensemble of given parameters (a, 6) as defined by ( ll.4p . Let us 
emphasize the fact that several different ensembles correspond to one set of parameters 
a, 6 in fOj) . 

The zoo of the known Meixner ensembles includes obviously the Gaussian ensembles 
which are described by their cumulant transform C\\x6 + y(tr#) 2 + ytr(# 2 ) with the 
familiar GOE, GUE and GSE (see j2] or [31]) corresponding to c 3 = 1 and C\ = c 2 = 
among them. The Gaussian ensembles yield a family of Meixner ensembles corresponding 
to a = b = in (11.41) . but the existence of non Gaussian ensembles for a = b = and n > 3 
has not been proved or disproved. Another important example is a Bernoulli ensemble: 
if /i m is the distribution of the projection on a uniformly random m dimensional space of 
the Euclidean space of dimension n, a Bernoulli ensemble is any mixing of Ho, . . . , /i n ; for 
n = 1 we get back the ordinary Bernoulli distribution. The Bernoulli ensembles can be 
seen as the only distributions of random projections P such that UPU* ~ P for U G K, n ^. 
The convolution of N identical Bernoulli ensembles provides the binomial ensemble, and 
randomizing iV by Poisson or negative binomial distribution yields the Poisson and the 
negative binomial ensembles. The construction of the full family of Meixner ensembles 
with 6 2 = 4a or with b 2 < 4a is available only for n = 2 and their extension to n > 3 is a 
challenge. 

The paper is organized as follows. In Section [2] we describe Laplace transforms of 
Meixner ensembles. In Section E] we give examples of Meixner ensembles. In Section H] we 
derive the system of PDEs that determine the Laplace transforms of the general n x n 
Meixner ensembles. The general solution of this system is elusive for n > 3. In Section 
E] we use the system of PDEs to show that under a natural integrability assumption the 
examples of Meixner ensembles from Section [3] exhaust all 2 x 2 Meixner ensembles. In 
Section |6] we collected additional material: Proposition 16.11 answers in negative a ques- 
tion raised in [HI Remark 5.8] for 2x2 random matrices with finite Laplace transform. 
Proposition 16.31 gives a series expansion for a Laplace transform of a random projection 
which can be used to derive (some) solutions to the system of PDEs. Section 16.41 makes 
some comments about the Jordan algebra context of the problem. 
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2. Laplace transforms 

In this section we rewrite Meixner property (11. 2\\ in terms of the Laplace transform. We 
will work simultaneously with symmetric, unitary and symplectic ensembles, thus with 
probability laws on H njj g which are invariant under the action X y UXU* where U is 
in the group K n $. We call ft = 1,2,4 the Peirce constant (see [16], page 97], where ft is 
denoted as d and n is denoted by r). The space H U)/ a is a real vector space equipped with 
the Euclidean structure defined by the real inner product 

(x,y) i y (x\y) := U(tr(xy)) 

(the real part of the trace is needed for the symplectic case ft — 4 only). For an H n ^- 
valued random variable X, the Laplace transform is a mapping from H n>/ 3 into (0,oo], 
defined by 

HO) = E(exp((0|X»). 

Let 

(2.1) 6x = int{6 : L(0) < oo}. 

(We will also use notation M when fi is the law of random variable X.) Throughout this 
paper we consider only random matrices X with values in H nj/ g with non-empty Ox- 

We need a result of linear algebra that we are not going to prove. Consider the space 
L s (M n ^) of symmetric endomorphisms of the Euclidean space H n>i g. To each y £ H„ 5j g we 
associate the elements y ®y and F y of L s (H„ ij g) defined respectively by 

ht->(y®y)(h)=ytT(yh), h^F y (h)=yhy. 

They provide important examples of elements of L s (H nj/ g). Now, L s (H n ^) is itself a linear 
space, and the result that we are going to admit as a black box is the following (see [HI 
Lemma 6.3] and [271 Proposition 3.1] for a proof): 

Proposition A. There exists a unique endomorphism \1/ of L s (H ni/ g) such that for all 
y £ H ni/ g one has ty(y ® y) = F y . Furthermore 

(2.2) *(P«) = fy®y+(i-0P»- 

With this notation we have the simple but crucial result. (Note that this result could 
remove the hypothesis of invariance under rotations.) 

Proposition 2.1. Let fi be an ensemble on H n / j such that its Laplace transform L is 
finite on an open non empty set. The following are equivalent: 

(i) The ensemble \x is Meixner with real parameters A,B,C, i.e. (1 1.2ft holds for 
X ~ fj,. 

(ii) With k(d) = In L(9), we have 

(2.3) 2(1 - A)^(k"(6))(l n ) = AA(k\6)f + 2Bk\6) + Cl n 
for all Be G x - 
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Proof. For showing (i)=^(ii) we first observe that 

(2.4) E ((X - Y) ® (X - Y)e <9|s) ) = 2L"(6)L(6) - 21/ (9) <g> L'{6). 
To see this, note that since X, Y are i.i.d, the left hand side of (12. 4 j) is 

2L(#)E (X <g> Xe< e l x >) - 2E (Xe< 9 l x >) ® E (Ye^) . 
The same reasoning gives 

(2.5) E (S <g> Se (e|s) ) = 2L"(6)L(6) + 2L'(0) <g> L'(0). 

Since E (Se< l s >) = 2L(9)L'(9), applying * to flU) and CTl . from (O) we get 

(2.6) 2(1 - AML"(9))(l n )L(9) = 2(1 + A){L'{6)f + 2BL(9)L'{9) + C(L(0)) 2 I n . 

Since L' = fc'L and L" = (k" + k' g> fc')L, this implies fl2T3jl . 

The proof of the converse implication is omitted. (Compare [231 Section 1.1.3].) □ 

In the above proposition, the use of the endomorphism \I/ may seem surprising. For 
a random variable X valued in a linear space E, the link between E(X £g> X) and the 
Laplace transform L of X is easy through the second differential of L. However, if E is an 
algebra, like the space of square real matrices, relating E(X 2 ) to L is difficult. Here our 
reasoning was based on Px(^4) = XAX and thus Px(In) — X 2 . Some other facts about 
\I/ are known; it is a symmetric operator on L s (H n)J g) with only two eigenspaces: the one 
generated by all the <g> y + ¥ y for the eigenvalue 1 and the one generated by all the 
|y <8> y + Pj, for the eigenvalue —/3/2. Details are in [29J Section 5] where the dimensions 
of these two subspaces are also computed. 

Remark 2.2. Note that if A = 1 then X is degenerate. Indeed, from (12. 3p we see that 
A — 1 implies k'(9) = const. 

We now show that Meixner ensembles are preserved under convolution power. 

Definition 2.1. Let /x be a probability measure on a finite dimensional linear space such 
that its Laplace transform is finite on a non empty open set. The J0rgensen set A(/i) of 
/i is the set of a > such that there exists a probability measure \i a with L^ a = L° and 
@ Mq = ®m> see e -g- P31 Page 767]. 

For instance A(/i) = (0, oo) if \i is infinitely divisible and A(/x) is the set of positive 
integers if /x is the Bernoulli distribution on {0, 1}. 

Proposition 2.3. Suppose 0^ ^ 0. If p, is a Meixner ensemble with parameters A ^ 
1,B,C and if a G A(/x), then fi a is a Meixner ensemble with parameters 

A a = A/ (A + a(l - A)), B a = aB/(A + a(l - A)), C a = a 2 C/(A + a(l - A)). 
Proof. Since k a (9) := k^O) = ak^{9), multiplying (j27J]) by a(l - A a )/(1 - A) we get 

2(1 - A a MK(9))(I n ) = {k'Mf + 2B ^:^ a) k'M + Ca <f-£°\ 
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Solving the equation 

A(l - A a ) 
a a(l - A) ' 

we get the formulas for A a , B a , C a such that (I2.3P holds. □ 

It will be convenient to have a version of Proposition 12.11 for standardized ensembles. 

Corollary 2.4. Suppose E(X) = ; E(X 2 ) = I n . Denote k(6) = lnE(e^). // fOD 
holds, then 

(2.7) *(k"(e)){I n ) = I n + 2bk'(9) + 4a(A;'(#)) 2 . 

Conversely, if the logarithm of the Laplace transform of X satisfies (\2.7\i ( together with 
the initial conditions) then (11.4j) holds, and E(X) = ; E(X 2 ) = I rt . 

Proof. We apply Q with A = 2o/(l + 2a), 5 = 26/(1 + 2a), C = 2/(1 + 2a). 

□ 

3. Examples of Meixner Ensembles 

In this section we give examples of Meixner ensembles. In Theorem 15.11 we will show 
that the examples exhaust the family of all Meixner ensembles on with G M ^ 0. 

3.1. Bernoulli and binomial ensembles. The Bernoulli ensembles (which are different 
than the Bernoulli random matrices of independent two- valued entries) will be our basic 
building blocks for more complicated ensembles. Suppose P m is the orthogonal projection 
onto the random and uniformly distributed m-dimensional subspace, i.e. P m = UI mjn U* 
where I m ,n is the n x n matrix with m ones on the diagonal followed by n — m zeroes, 
and U is uniformly distributed on the group K n ^. 

Definition 3.1. Denote by /i m the distribution of P m . A Bernoulli ensemble with pa- 
rameters q±, . . . , q n > that add to at most one, is the law 

At = g <5o + <?i/ii H h g„_i/u n _i + q n 5i n , 

where q = 1 - {q x H \-q n ). 

In terms of random variables, Bernoulli ensemble can be realized by the random variable 



(3.1) X 



with probability q = 1 — (gi + • • • + q n ) 
P x with probability qi , 

I n with probability q n . 



Proposition 3.1. A Bernoulli ensemble is a Meixner ensemble with parameters A = —1, 
B = 2, C = 0. 
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Proof. We note that for any pair of projections 

(3.2) (p_Q)2 = 2 (p + Q)_(p + Q) 2 . 

Applying this algebraic identity to random projections X and Y, where X is given by 
(13. ip and Y is its independent copy, we see that (jl.2p holds. □ 

We note that for the Bernoulli ensemble formula (12. 6p takes a particularly simple form: 
its Laplace transform L B {6) = E(exp((0|X))) satisfies 

(3.3) *{L" B {9)){I n ) = L' B {9). 
The corresponding version of (12. 3p is 

(3.4) *(k B (8))(ln) = k' B (9)(I n - k' B {9)). 

The fact that all Bernoulli ensembles share the same A,B,C will be convenient for 
calculations. But it will be easier to classify the resulting Meixner laws into the familiar 
families from the standardized form that appears in formula (11.41) . To do so, we find the 
moments of X. We first note that by rotation invariance, E(Pfc) = cl n , so Etr (P&) = nc. 
But tr (P fe ) = tr (l k)n ) = k, so 

(3.5) E(P k ) = E(P 2 k ) = -I n . 

n 

Therefore, with q — \{qx + 2g 2 + • • • + nq n ), we have E(X 2 ) = E(X) = ql n . 

(1-3) V 

;i - 2q)S 



If q 0, 1, the standardized matrix is X = .} _ (X — ql n ), and we get 

V'2( 1 -'?) 



(3.6) E((X-Y) 2 |S)=4 I n + 



2^(1 - ?) 4 



So (ll.4p holds with parameters a = —1/4 and b = (1/2 — q)/ \J q(l — q) . 

An intrinsic characterization of Bernoulli ensembles with an elementary proof is as 
follows. 

Proposition 3.2. // X with values in H nj/ g has a law invariant under rotations and 
X = X 2 then its law is a Bernoulli ensemble. 

Proof. Since X is symmetric, it is an orthogonal projection. Denote by c?(X) the dimension 
of its image. Invariance under rotations means that for any bounded measurable function 
/, we have 

E(/(Z7XET) - /(X)) = 
for all U G JC n> p. If d — 0,1, ... ,n and if / is zero outside of matrices of trace d, then 

E(f(UXU*) - f(X)\d(X) = d) = E(f(UXU*) - f(X)) = 0. 

This implies that the conditional law C(X\d(X) = d) is invariant under rotations, so it is 
a law of Pd- Then X is a Bernoulli ensemble with parameters q m = Pr(d(X) = m). □ 

This implies the following converse of Proposition 13.11 
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Proposition 3.3. Suppose that a law \x belongs to a Meixner ensemble with parameters 
A = — 1, B = 2, C = 0, and that the first four moments are finite. Then /i is a Bernoulli 
ensemble. 

Proof. Let X, Y be independent random variables with law fi, and let S = X + Y. Since 
(X - Y) 2 = 2X 2 + 2Y 2 - S 2 , and since the law of (X, S) is the same as (Y, S) from (TP]) 

we get 

E(X 2 |S) = ^S. 

From E(X|S) = |S we see that 

(3.7) E(X 2 |S) = E(X|S). 
From this, we deduce that E(X) = E(X 2 ), and that 

E(X 2 S) = E(XS), 

which implies that E(X 3 ) = E(X). 
Using (13. 7p again, we have 

(3.8) trE(X 2 S 2 ) = trE(XS 2 ). 

By rotation invariance, the moments E(Y J ) = Mjl n are given by real numbers, so 
E(tr(X r Y s )) = nM r M s = tr (E(X r )E(Y s )). Therefore, implies that 

tr E(X 4 ) + 2 tr (EX 3 EX) + tr ((EX 2 ) 2 ) = tr E(X 3 ) + 2 tr (EX 2 ) 2 + tr (EX 2 EX). 

Since we already proved that M X = M 2 = M 3 , we see that nM 4 = tr E(X 4 ) = tr E(X 3 ) = 
nM l . 

Let Ai, A 2 , . . . A ra be the (random) eigenvalues of X. From the equality of the first four 
moments of X, we see that 

rt 

E(A j (l - Aj)) 2 = tr E(X 2 (I n - X) 2 ) = n(M 2 + M 4 - 2M 3 ) = 0. 

i=i 

Therefore, all eigenvalues of X are either or 1, and X 2 = X. So by Proposition 13.21 the 
law fi of X is a Bernoulli ensemble. □ 

The J0rgensen set of a Bernoulli ensemble is of interest due to connections with free 
probability which allows continuous values for the analog of parameter N. For n = 2, we 
will show that the J0rgensen set is IN. The following sufficient condition shows that this 
is a "generic case" . 

Proposition 3.4. If a G A(/i) and fi is a Bernoulli ensemble with parameters qi,...,q n 
on H nt p, then (qo + q\z + • • • + q n z n ) a is a polynomial in variable z. In particular, na is 
an integer. 

Proof. For 6 S = diag(s, s, . . . , s), the Laplace transform is L{6 S ) = J2r=o <?re rs = rij( eS — 
Zj) mj , where Zj G C are the distinct roots of the polynomial go + Qi z + • ■ ■ + q n z n taken 
with their multiplicities rrij. Since the Laplace transform must be an analytic function 
on its domain, we see that a must be rational, and that if a = p/q with relatively prime 
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p, q e IN, then q must be a common divisor of all multiplicities mi,m 2 , ... of the roots. 
So q divides also their sum m\ + m 2 + ■ ■ • = n. □ 

3.1.1. Binomial ensemble. The real, complex or quaternionic binomial ensembles are mea- 
sures on H n>i g which are parametrized by integer iV and a discrete probability law on 
{0, 1, . . . , n}. We will follow the tradition of not listing the probability of among the 
parameters. 

Fix integer iV and non- negative numbers qi,...,q n with Qi + • • • + Q n < 1- Let 
Xi, . . . , Xjv be independent random matrices with the same Bernoulli distribution (13. ip . 

Definition 3.2. The binomial ensemble Bin(iV, qi, . . . , q n ) with parameters N = 1, 2, . . . 
and qi, . . . , q n , is the law of X = J2j=i ^-j- 

Proposition 3.5. A binomial ensemble with parameter N is a Meixner ensemble with 
parameters A = -1/(2N - 1), B = 2N/(2N - I), C = 0. 

Proof. This is a special case of Proposition ^. 3[ applied to the law \i of Xi with parameters 
described in Proposition 13.11 and to a = N. 

□ 

For standardization, we will need to know that 
(3.9) E(X) =iVgand Var(X) = Nq(l-q). 

From f TL5|) we then get that (H3D holds with a = -1/(4N) and b = (l/2-q)/y/Nq(l - q). 

3.2. Poisson ensemble. Our Poisson ensembles have parameters Ai, A 2 , . . . , A n > with 
A = Ai + ■ • ■ + A n > 0. 

Definition 3.3. Let be a Poisson real random variable, Pr(A^ = j) = e~ x X> j = 
0, 1, . . . , and let q\ = Ai/A, . . . , q n = A n /A > 0. Let X 1; X 2 , ... be independent Bernoulli 
matrices with the same parameters q±, . . . , q n , and Xo = 0. The Poisson ensemble is the 
law of X = Y.k=o^k- 

We will use notation Poiss(Ai, . . . , A ra ). 

Proposition 3.6. The Poisson ensemble is a Meixner ensemble, as it satisfies (ll.2p with 
parameters A = 0, B = 1, C = 0. 

Proof. Let Lb{0) be the Laplace transform of the corresponding Bernoulli ensemble. Then 
the log of the Laplace transform for the Poisson ensemble is k(6) = \{Lb{0) — 1). From 
(13. 3 p we get 

*(k"(6))(I n ) = k'(8), 

which is (TO]) with A = 0, B = 1, C = 0. Therefore §£2$ holds by Proposition El □ 
For standardization, we will need to know that, with A = (Ai + 2A 2 + • • • + n\ n )/n, 



(3.10) E(X) = Al n and Var(X) = AI n . 

Indeed, E(X) = E(E(X|JV)) and 

Var(X) = E(Var(X|N)) + Var(E(X|N)) = E(Nq(l - q))I n + Var(Nq)I 
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One can then check that (11.41) holds with a = and b = It is also easy to see 

that if the law of X is Poisson with parameters Ai, . . . , X n > 0, then X = Ylm=i -^m 
is the sum of independent random variables X m from Poisson ensembles with parame- 
ters (0, . . . , 0, A m , . . . , 0). Furthermore, real random variable tr (X) has the compound 
Poisson law with the law of summands given as X^=o ^x$k- 

We remark that the Poisson model with exactly one parameter Ai ^ appears explicitly 
in [121 P a g e 638], as part of the construction of matrix models for all free-infinitely divisible 
laws; see also [7]. However, regression properties of the model were not analyzed. 

3.3. Negative binomial ensemble. We now use Bernoulli ensembles to construct the 
negative binomial ensemble. Let N be a negative binomial random variable with param- 
eters r > 0, < p < 1, i.e. 

(3.11) P(N = j) = ^t^ PV, J = 0, 1, 2, • • • , q = 1 - p. 

Definition 3.4. The negative binomial ensemble NB(r, g 1; . . . , q n ) with parameters r > 
and qi, . . . , q n > such that qi + - — h q n < 1, is the law of the random sum 

TV 

x = ^ x fe , 

k=0 

where N has distribution (13.1 ip with p = 1 — (q x H — • + q n ), X 1; X 2 , . . . , are independent 
Bernoulli ensembles with parameters qi/(l — p), . . . ,q n /(l — p), and X = 0. 

Proposition 3.7. The negative binomial ensemble is a Meixner ensemble, as it satisfies 
(ll.2j) with parameters A = ^rj, B = ^+i> (7 = 0. 

Proof. Let L B {0) be the Laplace transform of the corresponding Bernoulli ensemble. The 
generating function ~E(z N ) = p r /(l — qz) r gives 

(3.12) L(9) = E(e< e l x >) = p r (1 - qL B {6)y r , 
so 

l-qL B {6)=p{L{6))- l /\ 

Differentiating, we get 

qL'JB) = ^L'(6)(L(6))^ 1+r ^ r 
r 

and 

qL" B {9) = ^L"(6)(L(6))- {1+r V r - g£i±rl i/(e) <g> L'{B){L{B))-^ + ^' r . 
From (I3.3p . we get 



We re-write this as 



2r , T . /W///1WT 2(r + l) /w/ _ 2 . 2r 



*(L"(#)(I n ))L(£) = -^_^(L'(#)) 2 + -—L'(9)L(6), 



2r + l v v /v " y ' 2r + l v K " 2r + 1 
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which is ( 12. 6ft with A = t^tj, -£> = 2^+1 ' C = 0. Therefore (II. 2p holds by Proposition 

El r r ^ 

For standardization, we will need to know that with q = {qi + 2g 2 + • • • + nq n )/n, 

(3.13) E(X) = rh n and Var(X) = ^| (p + q) I n . 

p p 2 

Indeed, E(X) = E(E(X|JV)) = q/(l - p)E(iV)I„ = qr/pl n and 
Var(X) = E(Var(X|N)) + Var(E(X|N)) = E(Nq(l - q))I n + Var(Nq)I n 



q(l--^)-I n + q 



-2 



1 — p p n (1 — p)p 2 n 
Then from (11.51) we see that (1 1.4ft holds with a = 4- and b = — £=£=. 

2- v /rg(p+g) 

3.4. Gaussian ensemble. Since X — Y and S are independent for any Gaussian inde- 
pendent identically distributed pair X, Y, formula f ll.2j) holds with A = B = and C — 2. 
The requirement of rotational invariance reduces the choices of the Gaussian law to the 
three-parameter family: 

(3.14) fc(0) = ci tr + c 2 ( tr #) 2 /2 + c 3 tr (0 2 )/2 

with C3 > and nc2 + C3 > 0. To see this, without loss of generality we take c% = 0. Let 
61, . . . , 8 n be the eigenvalues of 6, and consider the matrix of the quadratic form 02(^1 + 
• • • + n ) 2 + cz{6\ H — • + 8^) whose characteristic polynomial is (z — (nc 2 + c 3 ))(z — c 3 )" _1 . 
This shows that c 3 > and nc 2 + c 3 > 0. 

The explicit construction is to start with auxiliary GUE/GOE/GSE matrix Z and 
independent real standard normal (. 

Definition 3.5. For c\ G R, c 3 > and c 2 G E such that nc 2 + c 3 > 0, the Gaussian 
Meixner ensemble is the law of 



X = Z--tr(Z)I n J + Wc 2 + -CIn + ciI n . 



n / v n 



Proposition 3.8. The logarithm of the Laplace transform of the Gaussian Meixner en- 
semble is given by ( 13 . 14j) . 

Proof. Since (0|Z-atrZ) = (0-atr0|Z> and tr ((0 - ± tr 0) 2 ) = tr 2 - i( tr 0) 2 , we see 
that the answer follows from the well known GUE/GOE/GSE formula Eexp(tr (0Z)) = 
exp(tr0 2 /2). □ 

3.5. Gamma ensemble. The remaining types of Meixner ensembles will be constructed 
only for matrices of size n = 2. One difficulty we encounter is lack of continuity, which 
we now explain. 

It is natural to expect that gamma ensemble arises as lim p _ 5 .o£'Xp of a sequence of 
negative binomial ensembles with varying parameter p while r and the ratios qj/(l — p) 
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are kept fixed. However from (13.121) we see that 

lim Ee <e|pXp) = \wvp T (l- qL B (6)Y r = (l- tr (BL' B (0))- r = (1 - gtr 0)) _r . 

So this is a trivial ensemble of the form £I n with real gamma-distributed £. 

We now show that there are non-trivial gamma ensembles for n = 2. This construction 
is based on a more detailed analysis of the system of PDEs that arises from (12. 7p . 

Definition 3.6. The Gamma ensemble on H2./3 with parameters p > (3/2, c > lis defined 
by its Laplace transform 

(3.15) E(e (e|x) ) = (l - W 1 + (3 tr 9 + (3 tr 2 6 + 4 det 6 s ) 

defined on 

X = {6 e H 2)/9 : 1 -4c v / l + /3tr# + /3tr 2 # + 4det# > 0}. 

Of course, this definition requires a proof that the required law on H 2)/ g exists. The 
Gamma ensemble on H2, J g is constructed by choosing the appropriate law on K/ 3 " 1 " 2 , and 
arranging the corresponding real random variables into the random matrix. The construc- 
tion is based on the laws analyzed by Letac and Wesolowski [28l Theorem 3.1]. According 
to this result, for p > (3/2 and c > 1 there is a probability measure u P;C (dx) on the open 

Lorentz cone fig = {x G R^ 2 : x > JYlj=x X V\ with the Laplace transform 



(3.16) 



0+1 



3 5Z s j x 3 



i> PtC (dx) 



1 - 2cs + si - 



i=i 



Proposition 3.9. Let (£ , ■ ■ ■ , 6/3+1) have joint distribution z/ PiC 
wii/i parameters p, c is: 



The gamma ensemble 



(i) /or /3 



(3.17) 



(ii) /or = 2, 



(3.18) 



X 



X 



v/rr^o + ei 

6 



6 



6 - *6 



6 + if 3 



where i = 
Ciii) /or /9 = 4, 



(3.19) 



X 



6 + ^3 + ^4 + ^6 



v/T+Tfo - 6 

where i,j,k are the standard quaternion basis. 
Proof. We give the proof for (3 = 4. Writing 

«o + s i s 2 + is 3 + js 4 + ks 5 



(3.20) 



e 



S2 — isz — S4 — ks 5 



so - Si 
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we have (0|X) = 2 v / 5s fo + 2 £j =1 s&, so (l3~T5l) follows from (KTIBjl . □ 

Proposition 3.10. ^4 gamma ensemble on H 2/ g ifii/j parameters p > (3/2, c > 1 is a 
Meixner ensemble with parameters 

A = B = 0, C = 0. 

1 + 2p' 

T/ie mean and tae variance are E(X) = 2pcy/l + /3 12, Var(X) = 4pc 2 (l + /?)I 2 . For the 
standardized version, (11 ,4ft /ioWs wrat/i a = 1/ (4p) and 6 = 1/ y/p. 

The proof relies on the system of PDEs derived in Theorem l4.lt it appears after Propo- 
sition 15.61 

Remark 3.11. Somewhat more generally, one can define gamma ensembles on H n / g as 
the ensembles with Laplace transform 

(3.21) L{9) = (1 - ctr0 + (/3(n - 1) + 2)(trfl) 2 - 2ti6 2 )- p . 

Such an ensemble is a Meixner ensemble, as it follows from Theorem 14.11 that (ll.4p holds 
with a = l/(4p), b = 1/y/p. Since this is only a one-parameter subset of the possible 
solutions, we do not pursue this construction further. 

3.6. Hyperbolic ensemble. Recall that the Bessel functions and the modified Bessel 
functions are [381 $3.1 (8), $3.7 (2)] 

(3 , 2 ) ^)=E ( ;™:;: . u^±^f-. 

*— i n\ 1 in + v + 1) ^n!Im + i/ + l 

n=0 v ' n=0 v ' 

Following [23], we work with differently normalized Bessel functions 

(3.23) J v {x) := T{v + l){2/x) v J v {x), 

(3.24) X v {x) := T{u + l){2/xY I„(x). 

which are entire functions of z e C for u > —1. 

Definition 3.7. Fix a > 0, < A < 1, p 6 R. Define G (— 7r,7r) by sin<^ = 
p/a/(1 - A) 2 + p 2 . Let ^ be the first positive zero of J^-i)/ 2 {x) i.e., t x « 2.40483, 
£ 2 = 7T « 3.14159, and £ 4 » 4.49341. 
Let 

6x = {0 G H 2>/3 : |6>! + # 2 + y?| < vr/2, ^ - 2 | < ^} , 

where 9i, 9 2 are the eigenvalues of 9. 

The hyperbolic ensemble on H 2)/ g with parameters a, A, p is defined by its Laplace 
transform: 

(3.25) E(e (e|x> ) = ((1 - A) cos tr 9 - p sin tr 9 + XJ^-iy 2 (V tr 2 # - 4 det 9 
on the set 0x- (When p = we also allow A = 1.) 



14 
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These ensembles are constructed again by finding an appropriate law on R /3+2 and 
arranging the real random variables into a random matrix in H 2ij g. We begin with an 
auxiliary fact about Laplace transforms of probability measures on product of Euclidean 
spaces. 

Proposition 3.12. Let E and F two Euclidean spaces, and S C E, T C F be open 
non-empty sets. If f : S — >■ R g : T — >■ R are such that 



are Laplace transforms of probability measures for some a > and also for a — 1, then 
for < A < 1, and (s, t) G S x T, function 

1 



is a Laplace transform of a probability measure H a (dx, dy) on E x F. 

Proof We denote by j a the gamma distribution on (0, oo) with shape parameter 1 and 
scale parameter a > 0. Fix < A < 1, and let ji a and u a be two probabilities on E and 
F such that the integrals 



are finite on the maximal open sets S C E and T C F respectively. (Denote \i = \i\ and 
v = z/i.) 

Consider the set of independent variables U, A, B, Xi, X 2 , . . ., Yi,Y 2 , . . ., Nx, Ny 
where U ~ 7 Q , A ~ fi a , B ~ u a , Xj ~ /x, Yj ~ u, and t >->■ N x (t) and t h-> iVy(t) 
are two standard Poisson processes. Define the random variables X and F on E and F 
respectively by 



(3.26) 



1 1 



f"(sY g°(t) 



((l-\)f(s) + \g(t)T 




N X {XU) N Y ((1-\)U) 



X = A + E x ^ r = 5 + E ^- 




We have 



E ( e (s\x) + (t\Y)\u} 



1 



e-^/R-^e-^-^^iT*)- 1 ). 



f a (s)g«(t) 



Thus 




1 





□ 
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We now apply this result to the following Laplace transform, with the notation of 
Definition 13.71 

Proposition 3.13. For every a > 0, m > 1, /? > —1, p G R ; and < A < 1 there exist 
a unique probability measure H = H a ^ t \^ m onEx R m with the Laplace transform 

(3.27) 

/ / exp(sx + (t,y))H(dx, dy) = -, 

J iiRxR m ' ((1 - A)(cos(s) - psin(s)) + XJ^-i)/ 2 (\\t 



for s E R, t E R m such that \s + <f>\ < n/2 and \\t\\ < 1$. 

Proof. According to [381 $15.27], Bessel function J{B-\\n has simple real zeros that come 
in opposite pairs. Arranging the positive zeros in increasing order, < j i < j% < . . . , 
from [381 $15-41 (3)] we get 



J i e- 1)/2 (z) = l[(i-z 2 /j 



k=l 



so for < r < ji 



Hi 



oo ^ 



is a Laplace transform (of the convergent series of independent symmetrized Gamma 
random variables). Therefore, by Schoenberg's theorem [361 Theorem 2], for every m 
there exists a measure v a on Borel sets of R m such that for ||i|| < £g, 

1 



(J, 



(P-l)/2\ 



\t\ 



exp(t,y)u a (dy). 



Similarly, since cosi 



1 + P 2 ) 1 ^ 2 , we have 



cos(s) — psm s 



cosfs + 



cos « 



^ l-4(s + 0) 2 /(7r(2fc-l)) 2 
11 1 -40 2 /(tt(2A;- l)) 2 ' 



k=l 



so there is a probability measure \i a such that for \s + <p\ < tt/2, 

^ e sx p a (dx). 



[cos s — p sin sj - j R 

Proposition 13. 121 shows that there is a probability measure H a ^^ P)m on R x R m such that 
( l3~27j) holds. □ 

Proposition 3.14. Let (£ , • • • , ip+i) have joint distribution -ff a ,/3,A,p/(i-A),/3+i- The hy- 
perbolic ensemble on H 2j /3 with parameters a, A, p is: 

(i) for f3 = l, 

6 ^o-6 



16 
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(ii) for = 2, 



X 



6 - i& £0 - 6 



where i = 
(iii) for (3 = 4, 



1, 



(3.28) 



X 



£0 + 6 6 + ^3+^4 + ^5' 

6 - ^3 - jU - K$ £0 - 6 



where k is the standard basis of quaternions. 



Proof. To get the Laplace transform, we use (13.201) (or its appropriate modifications 
for < 4) with (0|X) = 2s £ + 2^5^. From (13371) with s = 2s and t = 
(2si,2s 2 , • • • ,25,3+1) we get 

(3.29) E(e^=° 2s ^) = ((1 - A)(cos(2 So ) - psin(2 So )) + AJ- (/3 _ 1)/2 (||t||))" a , 

which gives ( I3.25p . □ 



Proposition 3.15. The hyperbolic ensemble with parameters a > 0, < A < 1, p G R is 

a Meixner ensemble with parameters 



A 



5 = 0, C 



2a 2 {p 2 + 2) 



l + 2a 2a + 1 

The mean and the variance are E(X) = poU.^, Var(X) = a(l + p 2 )l2- For the standardized 
version, ( II ,4p holds with a = l/(4a) and b = pj J a{\ + p 2 ). 



We remark that the "exceptional" hyperbolic ensemble with parameters a > 0, A 



p = is also a Meixner ensemble with parameters A 



B = 0, C 



4a 2 



with 

mean E(X) = and with variance Var(X) = crf.2. The proof relies on the system of PDEs 
derived in Theorem I4.1( it appears after Proposition 15.61 



4. The system of PDEs 

In this section, we derive a system of PDEs for the Laplace transform of general Meixner 
ensembles. As previously, we consider simultaneously the real, complex, and quaternionic 
cases. 

We introduce some notation. For 9 G H n p and j = 0, 1, . . . , we consider (Tj(9) defined 

by 

n 

(4.1) det(I n + z6) = J2 a i( e ) zj - 

Recall that for = 4, the determinant if quaternionic matrices are defined only for 
hermitian matrices, see [TBI P a g e 29] and [IB] . Note that 

<7„ = 1, <7i(0)=tr0, a 2 {e)= l -{ti6) 2 - l -ti{6 2 ) 1 a n {6) = det 6 , 
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and that Oj = for j > n. We also adopt the convention that Oj = for j < 0. In 
general, <Tj(9) — &j{0\, 9 2 , ■ ■ ■ ) is the j-th elementary symmetric function of the eigenvalues 
61,62, ■ ■ ■ of G H nj/ g; this notation and most of our calculations do not depend on the 
dimension n. Denote by U the open subset of R n such that if (<7i, . . . , a n ) G U then the 
polynomial in x defined by ^™ =0 (— l) n ~ l (JiX n ~ l has only distinct real roots. 

The main result of this section is the system of PDEs that determines the Laplace 
transform of a Meixner ensemble on a nonempty open set 0x C H n ^. This system is 
written in terms of an auxiliary function g{(T\, . . . , u n ) which is defined on an appropriate 
non-empty open set Ux C U. The link between L and this auxiliary function g varies 
according to the fact that a or b are zero or not. To define the set Ux, we need 0x to be 
closed under conjugation {0 G ©x implies U0U* G ©x for all U G fC n p), and to have in 
its closure. Consider the n-dimensional real subspace D n C H„ >( g consisting of diagonal 
matrices, and let D Q n C D n denote the set of matrices with distinct eigenvalues. Since ©x 
is invariant under rotations, ©x H D n is a nonempty open set with in its closure. The 
set ©x H D® is obtained from ©x H D n by removing a finite number of hyperplanes, so it 
is also non-empty, and has in its closure. 

Denote by a : H„ j/3 -)■ R n the mapping ^ (o x {0), . . . , a n (0)). If G D° then by the 
Vandermonde's determinant, Be = {l n , 0, 2 , . . . , n ~ 1 } is a basis of D n . From formula 
(14. 10p below, we see that in that basis, the matrix representation of the derivative of o 
restricted to D n is triangular at G with <t = 1 on the diagonal, 

a Q ... 0" 

-ai{0) a ... 

<r 2 (0) -<ri(6) a ... 

(-l)"- 1 ^!^) (-I) n -V n _ 2 (0) (7 . 

Therefore, the Jacobian of a is 1 and o : Z^,° — > R™ is an open mapping. Thus, the set 

(4.2) u x = a(e x r\D°) 

is open, non-empty, and since c(0) = 0, it has G R n in its closure. 
To shorten the formulas, for i, j — 1, 2, . . . , n we write 

9i= and = a a • 
a<Tj d(Jid(jj 

Recall the convention that <jj = for j < or j > n. 



W{o)] Be = 
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For g : C/x - > R consider the vector 

9ll ~ E" s =2 (T r-+ S -25r S ~ (n ~ l)f </2 
+ 2^12 - Er, s =3 <Vh»-30« - (n - 2) §#3 



(4.3) 



EJ,s=l a r+s-l-j9rs — Er,s=j+1 Cr+s-l-jflVs ~ (" ~ j)^9j+l 



En— I _ /3 

r s= \ 0~r+s—n9rs ®n9nn ^dr, 

En 
r )S =l &r+s—n—l9rs 

Theorem 4.1. Let X be a random variable with values in H n .^ with the law invariant 
under rotations. Let L{9) = e k ^ e ' = E(e^' x ^) be its Laplace transform, and assume that 
©x 7^ 0. Suppose that E(X) = 0, E(X 2 ) = l n and that X satisfies (11.4)) with some a and 
b. Let £/x be the associated open subset o/R n , see (14. 2)) . 

(i) I/a^O, define g: ->• R by g(ax{9), . . . ,a n {9)) = exp{-4ak(6) - btr 6). Then 
g satisfies the system of PDEs 

r(b 2 -4a)g- 


K<7) 



Furthermore, 

(4.4) ^(tn, . . . , O ->• i, 



I cti=0 



-6 



as (o"i, . . . , cr n ) — > over Ux- 
(ii) Ifa = 0,b^0, define g : ^ R by g(<ri(0), . . . ,<r n (B)) = k(6) + ±tr6. Then g 
satisfies the system of PDEs 



Kg) = ib 



9i 

92 



Furthermore, 



(4.5) 



g(ai, . . .,a n ) -»> 0, 



dgjai, ...,a ri 



1 

^26 



as (<Ji, . . . , a n ) — > over £7x- 
(iii) If a = 0,b = 0, define g : C/x — ► R 5'(c r i(^), • • • , 0n(0)) = k(9). Then g satisfies 
the system of PDEs 

"1" 





o 
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Furthermore, 

(Aa\ ( \ , n dg(o- 1 ,...,a n ) 
(4.6) g(ai,...,a n ) -+ 0, — >0 

as (cti, . . . , a n ) —> over C/x- 

Conversely, let \x be an ensemble on such that M ^ 0. // one can find constants 
a,b that lead to one of the equations in (i), (ii) or (Hi), then the corresponding version of 
equation (j2.7p holds for all 9 £ G M; and hence fi is a Meixner ensemble with parameters 
a, 6. 

The left hand side O(g) of our systems of the PDE's resembles the system of PDE's for 
the Bessel function of matrix argument given in [211 Eqtn (18)] when /3 = 1, see also [UJ 
page 275 Eqtn (31)]. However, as we will see in Proposition I5.6[ in contrast to [21] our 
system has multiple analytic solutions. 



4.1. Proof of Theorem 14.11 The Cayley Hamilton theorem (for the case (3 = 4, see [T6| 
Proposition II. 2.1]) implies 

n 

(4.7) e n = J2(- 1 Y~ 1(T i{ d ) en ~ i - 

8=1 

We use also without proof the Newton formula 

m— 1 

(4.8) ma m {9) = ^(-lr^V^) tr (^" l ); 

we remark that under the convention that o m = for m > n, formula (14.81) holds true for 
all m > 1, i.e. also for m> n. Since trl n = n, it can be rewritten as 

m 

(4.9) (m - n)a m {6) = ^(-l)™- 1 "^^) tr (# m " 1 ). 

i=0 

We also need a differentiation formula. 
Proposition 4.2. Fza; n > 1. Taen /or a// m > 1, 

m— 1 

(4.10) o4(0) = ^(-l)— 1 " 4 ^^)^- 1 -. 

i=Q 

Proof. We prove (14.101) by induction with respect to m. Since o~[ = l n = o~q8°, the formula 
holds true for m = 1. Suppose ( 14.101) holds for some m > 1. With t,- = tr(0 J ), we 
differentiate Newton's formula (14.81) written for o~ m+ i and use the induction assumption. 
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We get 

(4.ii) ( m +i)Ci=H)>+ir 



E 1 

j=l s=0 



= (-l) m (m + l)# m + (m + 1) E(-l) m -V^ m - J ' + i2, 

where 

m m j—l 

r=- ;>j ir \irf,o" ! > + ^(-i) m ^wi-, j2(- i y~ 1 ~ S(T ^ j ~ 1 ~ s 

j=l j=l s=0 

m m s—1 

= - E(-l) m ^JCT^ m ^ + J2(-l) m ~ S m ~ S E(- 1 ) S-1 ~ i(7 iV J - • 
j=l s=l j=0 

Using Newton's formula again, we see that X^=o(~ l) s_1 ^ a jts~j = scr s , so R = and the 
formula follows from (14.111) . □ 

We now consider separately the following cases, in which we re-write (12. 7p by the 
indicated substitutions: 

(i) If a ^ 0, then f{9) = exp(-Aak(9) - bti (9)) solves 

(4.12) *(f"(9))(I n ) = (b 2 -4a)f(9)I n . 
The limits as 9 — > in 0x are: 

(4.13) /'(^^-feln- 

(ii) If a = 0, 6 ^ 0, then f{9) = k{9) + ^ tr (9) solves 

(4.14) M/(r(0))(I n ) = 26/'(^)I n . 
The limits as 9 — > in Bx are: 

(4-15) /(*)-► 0, /'(0)^I n . 

(hi) If a = 0, 6 = 0, then f(6) = k{9) solves 

(4.16) V(f"(9))(I n ) = I n . 
The limits as 9 — > in Gx are: 

(4.17) 0, f (0) 0. 

The remaining part of the proof consists of finding a suitable expression for /. 

Now the fact that the law of X is invariant under the rotations is equivalent to k being 
invariant under the substitution 9 i— > U9U* for all U G K, n> p, and is equivalent to saying 
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that in each of the three cases mentioned above, there exists a real analytic function g on 
Ux such that 

f{9)=g{a l {9),...,a n {9)) 

for all 9 G Ox with distinct eigenvalues. (This is Weyl's formula, see [3TJ Lemma 2.6.1].) 
Using (I4.10p we get 

n n m— 1 

(4.i8) f{9) = ^9^(0) = ^^(-lr- 1 ^)^- 1 - 4 - 

m=l m=l i=0 

The second derivative of / is trickier: we write f" = A + B with 

n n 

m=l i=l 
n 

m=l 

The calculation of ty(A)(l n ) is quick and gives 

n n 
m=l i=l 



For calculating ty(B)(I n ) we prove the surprising formula: 
Proposition 4.3. ^«J(I„) = f (m - 1 - n)a' m _ x . 
Proof. We split a'^(9) into a'^{9) =C m + D m where 





m— 1 
Cm = ( 


i-l 




i=0 


i=o 




m-2 


m-2-j 


(4.19) 


= £< 






i=o 


s=0 


(4.20) 


m— 1 

Dm = 


_l^-~l~j a ^ e m-l-jy_ 




i=o 





Using Proposition 1X1 we get 

/m-2-j \ 

(4.21) * J3 0" ® m - 2 -i- s (I n ) = (m — j — l)^" 2 , 
which leads to the calculation of \E , (C m )(I ri ) : 

m-2 

(4.22) *(C m )(In) = $3(-l) m -V,(m- j - l)9 m ~ j - 2 . 

j=0 
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Having in mind the calculation of ^(D m )(l n ) we write 



s-1 , s-1 



2 

i=0 i=0 



From Proposition we have ^(F x ) = (l - §) P x + fx <g> a; which leads to \^(P :r )(I n 
(1 - |)x 2 + fxtr (x) and to 



*(P I+S - P. - P,)ftO = (l - |) («y + yx) + g 'trM + gtrW 



Thus we can compute 



^([^ s ]')(in) = J f i - ^] Yye*-*- 1 + + j Y^e* tr (e-- 1 -*) + e 3 - 1 -' tr (^)] 

v 7 i=0 11 



i=0 i=0 
s-l 

3* / ns— X—i\ 



We apply this to s = m — 1 — j and compute \l/(.D m )(I n ) : 



(4.23) *(D m )(l n ) = (l - X>l) m ^(m - 1 -j)9 m -^ 



3=0 

m—2 m—2—j 



+ ^Y(-l) m ~ 1 ~ j cr j ^tr(^ m - 2 ^'- i ) 

j=0 i=0 
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Note that we replaced m — 1 by m — 2 in the first sum of the right hand side. We obtain 



^ 



2^ v '' ~ JK " L 

3=0 

m—2 m—2—i 



rm-2-j 



J tr i 

o j=0 

m—2 



= f ^(-lr-V^m-l-j)^- 2 ^' 

j'=o 

„ m—2 

(4.24) +| ^(-l)^(m - 2 - i - n)* m ..2-i 

„ m—2 

(4.25) = | ^(-1)^(^1 - 1 .,)()■" 2 ' 

3=0 
q m—2 

\r •(>"' 2 \j- n)*, 

3=0 

m-2 

= (3 ™ - - n J2{-l) m - j VjO m ~ 2 ~ j 

3=0 

m — 1 — n , 
(4-26) = o' m _ x . 

In this chain of equalities, the first is obtained by combining ( 14. 22ft and ( 14. 23ft and reversing 
the summation in i and j in the last sum; (14.241) is obtained by applying (14. 9 p (replacing 
m by m — 2 — z); (I4.25P is obtained by changing the variable j = m — 2 — i. The last 
equality (14.261) uses the equality ( I4.10p . Thus the Proposition 14.31 is proved. 

□ 

Gathering all terms, using the change of variable m = j — 1 in ^{B) and the fact that 
a' = 0, the expression *(/")(I n ) = ty(A + B)(I n ) becomes, 

n n p n—l 

w'xin) = + f Etf - n W • 

m=l i=l i=l 

Since eri(0) = I n , the three equations (14 . 1 2 [) . (14 . 14 j) and (14. 16|) respectively become 

n n p n—l 

(4.27) EE^>* + f Etf - n )^°i = - 4o )^)> 

m=l i=l j'=l 

n n „ n—l n 

(4.28) gmiv' m Vi + 2 E( J ~ n )9j+i<^j = 26 E 9m<r' m (8), 

m=l i=l jr'=l m=l 
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n n p n—1 

(4-29) ££<W^ + f £(j - n)flf i+l0 J = ^(0). 

m=l i=l j=l 

Noting that a[, a' 2 , . . . , o' n is a basis of the algebra generated by 9 with distinct eigen- 
values, the corresponding systems of PDEs arise by comparing the coefficients at o~j in 
(I4.27H4.29"]) . respectively. To end the proof, we only need to find the coefficients Pj(m,i) 
in the basis expansion 

n 

(4.30) a' m a\ = Y^PA^M- 

3=1 

This is accomplished by the following formula. 
Proposition 4.4. For all r,s > 1, 

r-l 

(4.31) a' r (ey s (6) = Y,Wv^i- ] ], 

3=0 

where 

[f,g] :=f{9)g'{9)-f{9)g{9). 

Proof. From (14.101) we get <Jj +l = aj8° — 9a j. Since the algebra generated by 9 is com- 
mutative, aj +1 a' n — aW n+1 = ((Tj9 — 9a'^)o' n — a'j(a n 9° — 9a' n ) simplifies to 

(4.32) a 'j+l a 'n - a j a 'n+l = [ a 3> a n\- 

We now prove the proposition by induction. We assume that (I4.3ip holds for some fixed 
r and all s > 1. This is trivially true for r = 1 since a[a' s — a' s = [a , a s ]. If (I4.3ip holds 
true for some r > 1 and all s > 1, then by f)4.32p and the induction assumption, 

r-l 

°V+l°s = G 'r G 's+l + Wr,cr s ] = ^2[aj,a r+s -j] + [<r r ,a s ], 

3=0 

which shows that (14.311) holds for r + 1 and all s > 1. □ 
Proof of Theorem\4-l\ Cases (i)-(iii) correspond to equations (I4.12H4.16]) . The left hand 



sides of the resulting three equations (I4.27H4.291 are then re- written using (I4.30p . From 
(I4.3ip we see that coefficients Pj(m, i) are linear with respect to ai, . . . , o n . Rewriting 
(rCTTI) as 

r+s— 1 r— 1 

^P d {r,s)a' j = £ (Tr+s-l-jOj -^(Tr+s-l-jVj, 
j j=s j=l 

we get the explicit formula 

fcr r+s _i_j if max{r, s} < j < r + s — 1 < n, 
—°~r+s-i-j if 1 < j < min{r, s} and j < r + s — 1 < n, 
otherwise. 
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With this notation equations (14.27114. 291) become the equations asserted in Theorem 14. 11 
To derive the corresponding boundary limits, we use (I4.18P and the respective limits 

dm, iisd, or g^D. 

It is also clear that these calculations can be reversed. From the equations in listed in 
(i)-(iii), we see that (I2.7p holds for all 6 G 6x with distinct eigenvalues. By uniqueness 
of analytic extension, (I2.7P extends to the whole open set 6x- □ 

Remark 4.5. Note that to have a solution of the system in each of the three cases is 
not a guarantee of existence of an ensemble: the function L(9) still have to be a Laplace 
transform of a probability measure. For example, for n = 2 there are solutions of the 
system for all a < 0, but according to Theorem 15. lf i) they are the Laplace transforms of 
a matrix ensemble only if — l/(4a) G IN. 

Remark 4.6. Particular Meixner ensemble can be found from the system of PDEs by 
restricting g to be a function of 01,02 only. In other worlds, we look for an ensemble 
such that the Laplace transform L(6) is a function of tr9, ti9 2 . Such strategy leads us 
to the gamma ensemble on H njj g mentioned in Remark 13.111 and to Definition 13.51 of the 
Gaussian ensemble. This strategy leads only to trivial ensembles for b 2 ^ 4a. 



5. Meixner ensembles on 2 x 2 matrices 

The goal of this section is to use Theorem 14.11 to prove that the constructions from 
Section [3] exhaust all Meixner ensembles on ~H.2,/3. 

Theorem 5.1. Up to affine transformations, the only Meixner ensembles on with 
the Laplace transform defined on a non-empty open subset of H^,^ are of the six types, 



depending on the values of a,b from (11.41) : Re-worded! 



(i) If fi is a Meixner ensemble on with a < and G M ^ 0, then — l/(4a) = 
iV G IN and fx is a binomial ensemble Bin(N, qi, q 2 ) from Definition 1 3. 1?l for some 
91,52- 

(ii) If fi is a Meixner ensemble on with a = 0, b > and B M ^ 0, then is a 
Poisson ensemble from Definition \3. 3\ for some A > 0. 

(iii) If // is a Meixner ensemble on H2 i( g with b 2 > 4a > and M ^ 0, then fi is one 
of the negative binomial ensembles from Definition \3.4 



(iv) If /j, is a Meixner ensemble on H 2i /3 with a = b = and 6^ ^ ; then /i is one of 
the Gaussian ensembles from Definition ^. 5\ 

(v) If n is a Meixner ensemble on with b 2 = 4a > and M ^ 0, then fi is one 
of the gamma ensembles from Definition \3.6[ 

(vi) If fi is a Meixner ensemble on H 2i /3 with b 2 < 4a and @ M ^ ; then fi is one of 
the hyperbolic Meixner ensembles from Definition \3. 7\ 

(The seventh case of non-random ensembles arises from degenerate affine transforma- 
tions.) 



5.1. Laplace transforms of binomial, negative binomial and Poisson ensembles. 
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Lemma 5.2. IfP is a random projection in Bin(l, 1,0), then with aj = <x,(0), we have 
(5.1) E(e< e l p >) = e^% p . l)/2 

where X v is a Bessel function (I3.24p . 




Proof. To construct P, we choose a direction 



x 2 



at random by taking two independent 



standard normal random variables as the components of the vector. Note that this means 
that Xi = Zi + iZ 2 in the case ft = 2 or X x = Z\ + iZ 2 + iZ 2 + kZ^ in the case ft — 4, 
where Z\, Z 2 , Z%, Z4 are the standard real- valued independent normal random variables. 
The matrix representation of P is 



1 



\Xi\ 2 X 2 Xi 
X\X 2 I X 2 1 2 



iXxP + lXal 2 

Due to invariance under rotations, without loss of generality we take diagonal 
'0i 0' 



6 2 



so that 



|P) = tr (0P) 



l^i| 2 + l^ 2 | 2 



(^IXxi 2 + e 2 \x 2 \ 2 ) = e 1 u + 2 (i - u). 



Since |Xl| 2 , \X 2 \ 2 are independent gamma distributed with shape parameter ft/2 and scale 
parameter 1/2, the distribution of U = |Jfi| 2 /(|X L | 2 + |X 2 | 2 ) is Beta(/3/2, ft/2). Denoting 
by R = I n — 2P the corresponding reflection, this gives 



E ( e (0|R>) = e fli+ft, E ( e -2(9|p>) 
_Ji + 0> m 



r 2 (/3/2) J 



-2(u0 1+ (i-u)e 2 ) u m-^i_ u )m-i du 



r 2 (/?/2) Jo 



Substituting u — (1 +v)/2 we get 



(5.2) 



Lr(0) 



2^ 1 T 2 (ft/2) j 1 

The integral is expressed in terms of the Bessel functions in [38j $3.71 (9)]. 
Remark 5.3. Denoting by 9\, 6 2 G R the eigenvalues of 0, 

f/o(0i-0 2 ) if /? = 1, 



□ 



E(e< 9 l R >) 



sinh(6>i-e 2 ) 
01-02 



cosh(gi-6> 2 ) _ sinh(ei-6>2) 



L-02) 2 



L-02) 3 



if /3 = 2, 



if /3 = 4. 
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Combining this with the Laplace transforms derived in the proofs of Propositions 13.51 
13.6} and 13 .7\ we get. 

Corollary 5.4. (i) //X has Bin(N,q 1 ,q 2 ) law, then its Laplace transform is well 
defined for all 9 e Hg, and is given by 

(5.3) E(e<*l x >) = e 1 "*" ^ e" CTl/2 + giX (/3 _ 1)/2 [j^ ~ ^ + 

where q = 1 — q\ — q 2 ■ 
(ii) I/X has Poiss(X\, X 2 ) law, then its Laplace transform is well defined for all 9 e 
Hg, and is given by 



a 2 



(5.4) E(e^) = exp ^e^Z^y, -± - a 2 j + X 2 e ai - X 
where A = Ai + \ 2 . 

(hi) //X has NB(r,qi,q 2 ) law, then its Laplace transform is well defined for 9 in a 
neighborhood of in H 2jj g, and is given by 

(5.5) E(e<«) = ^- ^ • 

h - Ue«/2T ( „_ 1)/2 (V^ - <T2j + g 2 e<" J J 

Remark 5.5. A calculation shows that for a given b 2 > 4a > 0, parameter p = l—q\—q 2 
in ( 15. 5 p can range over the set 

{ P E(0,l):p/(2-p)<K/b<p}, 

while g 2 = (bp — k)/k G [0, 1 — p] and gi = 1 — p — q 2 are uniquely determined by p. (Here 
K = Vb 2 - 4a, 6 > 0.) 

5.2. Solutions of the system of PDEs for n = 2. We consider in more detail the case 
n = 2 with 

(5.6) C/ = {K<7 2 ) : Aa 2 <a 2 }. 

We require that the solutions extend by continuity to o~\ = a 2 = 0. 
The system of PDEs simplifies to: 

(5.7) - o- 2 £f 22 = ^g 2 + (b 2 - Aa)g, 

(5.8) 2g 12 + a ig22 = 0. 
We will use as the initial conditions 

(5-9) s(0 ,0) = l,^U = -fc 

Denote k = \J\b 2 — 4a|. 
Proposition 5.6. Consider the system ( 15. 7p . ( 15. 8 p w't/i initial conditions (I5.9P . 
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(i) For b 2 = 4a > all solutions are 

g(a 1} a 2 ) = 1 - bo x + + 4<x 2 ), 

where C is an arbitrary constant. Accordingly, 

(5.io) m «p(-*.(»>/»> 



;i-fcT 1 (6l)+C(^<T 1 a (9)+4<T 2 (9)))V»' ' 
(ii) For 6 2 > 4a a// solutions are 

(5.11) aK (7 2 ) = C ie KCTl + C 2 e- KCT1 + C 3 l { ^ 1)/2 (K^/a 2 /4 - a 2 ), 

where Tm_\y 2 is defined in (13.241) . and C\, C 2 , C3 are arbitrary real numbers such 
that C\ + C 2 + C3 = 1 and C 2 — Ci = 6/k. Accordingly, 



(5.12) 



L(0) = e -^/(4«) ^C ie KCTl + C 2 e- Kffl + C 3 %_ 1)/2 (nyJaf-Aa^j 



-l/(4a) 



which can also be written as 

l/(4o) 



(5.13) 

L(0) = e~ 6<Tl/(4a) ^(1 - A) cosh(K(7a) - ^sinh(/€<ri) + AZ (/3 _ 1)/2 (^-4, 2 )) 



where X & is an arbitrary constant. 
(iii) For 6 2 < 4a all solutions are 



g(o~i, cr 2 ) = (1 - A) cos Ka% - — sin ko\ + \J(p-i)/ 2 (kiJ o\ - 4a 2 ^ , 
where A an arbitrary real constant. Accordingly, 

,, , / 5 / \ -l/(4a) 

(5.14) L(0) = e" 6,Jl/(4a) ( (1 - A) cosKfn - - sin/coi + \J^ 1)/2 (kJ a 2 - 4<r 2 ) 



Proof. The solutions that depend only on <j\ satisfy equation g n = (6 2 — 4a)g with two 
linearly independent solutions: 

(5.15) e Kai , e~ Kai if b 2 - 4a > 0, 

(5.16) cos(kcti), sin (kctx) if 6 2 — 4a < 0, 

(5.17) 1, 0-1 if 6 2 -4a = 0. 

To find the solutions that depend also on o~ 2 , we first note that (15. 8 p implies that 

(5.18) 2g 1 + a 1 g 2 = C(a 1 ), 
so 

(5.19) 2g 11 + g 2 + a 1 g 12 = C'(a 1 ). 
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We now eliminate gn and g± 2 from the equations. Subtracting from ( 15. 19ft the linear 
combination of 2 times equation (15. 7p and 0i/2 times equation (15.81) we get 

(5.20) K 2 /4 - a 2 )g 22 - ^g 2 - (b 2 - Aa)g = -C"(<7i)A 

For fixed a%, we consider ( 15.20j) as a differential equation with respect to a 2 in the interval 
o 2 < cr 2 /4, since we need g in this domain, see (I5.6p . For this reason, we denote t = 
crf/4 — a 2 and, assuming b 2 ^ 4a, we denote 

ts \ C"(«7i) 



2(o 2 - 4a) 

Thus g 2 = —y', g 22 = y" and (I5.20p becomes 

(5.21) 2ty" + (l + /3)y'±2K 2 y = 0, 

where the sign is chosen according to the sign of 4a — b 2 . Substitution y(t) = v^~^'' A u(x) 
with x = 2«Vt converts (I5.2ip into the Bessel equation 




Using [Tj 9.6.18] or [TJ 9.1.20], after some calculation one verifies that all solutions of ( I5.2ip 
which are bounded in a neighborhood of are proportional to X^_iy 2 {K\/t) when b 2 > 4a 
and to J(p-i)/ 2 (n\fi) when b 2 < 4a, see (13.241) and ( I3.23p . Substituting g(a\,a 2 ) = 

K{a 1 )X^_ l)/2 {^al/A - a 2 ) back into (JEl we get that ([EHD holds if K' = 0. So we 
deduce that K is constant, and a direct verification shows that g solves also (15. 7p . (We 
repeat the same reasoning with J{p-i)/ 2 instead of X^iy 2 if b 2 < 4a.) 

Similarly, one works out the solutions for b 2 = 4a > 0. In this case, after taking 
y{t) = g(cr 1 , cr 2 ) + 2<j 2 C'(o\)I (1 + f3) we get equation ( I5.2ip with n = 0. After some work 
this leads elementary solutions proportional to g{a\,a 2 ) = — I" °i anc l ^° additional 
solutions g(ai,a 2 ) = (erf — Aa 2 f l /3 ^ 2 for (3 > 1 or g(a\,a 2 ) = log(a 2 — 4cr 2 ) for — 1, 
which are unbounded at the origin. 

To conclude the proof of (0), we note that the general solution that is defined at o\ = 
a 2 = is 

g(a h a 2 ) = C x + C 2 o x + C 3 (/3 o x 2 + Aa 2 ). 

The initial conditions determine C\ = 1, C 2 = —6, which gives (I5.10p . 

To conclude the proof of (Jn]), we note that the general solution bounded at the origin is 
(15. lip , and the initial conditions are satisfied when C\ + C 2 + C 3 = 1 and C 2 — C\ = bj n. 
This gives (l5TT2|) . 

The initial conditions similarly imply ( I5.14p . 

□ 



We are now ready to prove that gamma and hyperbolic Meixner ensembles are indeed 
Meixner. 
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Proof of Proposition VJ.IOX To calculate the moments, we differentiate (13.1 6p to get E(£ ) = 
2cp, E(£i) = 0, Var(£ ) = 2p(2c 2 - 1), E(ff) = 2p, and Efofo) = Ofor ^ ^ £ 2 . 

For the standardized gamma ensemble X with E(X) = and E(X 2 ) = I, from (13.151) 

we get 

(5.22) E(e <9|S:> ) = e-^ tre (l - tr 6/^ + D(f3 tr 6 2 + 4 det 6))~ p , 

where .D = l/(4pc 2 (l+/3)). This is (I5.10p with p = 1/b 2 , so property ( 11.2ft follows from the 
converse part of Theorem 14.11 For the standardized version, (ll.4p holds with a = 1/ (4p) 
and b = \j^fp. Formula for the parameters A,B,C is now re-calculated from formula 

(Q. 

Note that the admissible ranges of parameters are p > (3/2 and < D < l/(4p(l + 

□ 

Proof of Proposition \3.15l We compute the moments of X entrywise: differentiating (I3.29P 
we get E£ = pa, E£ 2 = o(l - A + p 2 (l + a)) and for j > 1, Ef 3 - = 0, E£ 2 = aA/(l + p), 
E^Cffl) = for l x ± £ 2 . 

For the standardized hyperbolic Meixner ensemble X with E(X) = and E(X 2 ) = I, 
from (13.251) we get 



Ee <e|x) = e - p ^ a/{1+p2) (1 - A) cos( tr 0/y/a(l + P 2 )) ~ P sin( tr 6/^a(l + p 2 )) 



+ A Jo9-i)/ 2 (V tr 6 2 - 4 det 0/ y/a(l + p 2 )) 



Substituting p = b/n = b/y/4a — b 2 , and 1/a = 4a, we get (I5.14p . So from the converse 



part of Theorem 14.11 we see that (11.41) holds with a = l/(4a) and b = pj ya(l + p 2 ). 
Formula for the parameters A, B, C is now re-calculated from formula (11.51) . 

□ 

5.2.1. Poisson case: a = 0, b ^ 0. The system of PDEs simplifies to: 

(5.23) 0U-O2022 = ^92 + 2bg 1 , 

(5.24) 2# 12 + ct 1 # 22 = 2bg 2 . 
We seek solutions such that 

(5-25) g(0,0) = 0, dgi ° 1,0) \ = \- 

v J yv > ; > ^ | CT1=0 2fe 

Proposition 5.7 (Poisson ensemble). Consider the system (15.231) . ( 15.24ft wi/i initial 
condition (I5.25p . TTie solution is 

<?K<7 2 ) = ^ ((1 - C)e 2bCTl + (2C - l) e ^X ( ^ 1)/2 ^6 V / a 2 -4^ - , 

where C is an arbitrary constant. Thus 
(5.26) 



L(0) = exp (-^ + ^2 ((1 - C)^ CT1 + (2C - l)e feCT1 Z (/3 _ 1)/2 (by/af - ^ 



C 
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Proof. The solution that depends on o~\ only is in [25J; in our notation, 

g(a 1 ) = C 1 + C 2 e 2b ^. 

To find the solutions that depend on both variables, we first note that f !5.24j) implies 

(5.27) 2g 1 + a 1 g 2 -2ag = C(a 1 ), 
which we differentiate with respect to <j\ to get 

(5.28) 2g n + g 2 + a x g 12 - 2ag x = C"(<ti)- 

We now eliminate gu, g\ 2 and g\ from the equations. To do so, from equation (15. 28ft we 
subtract 2 times equation (I5.23p . b times equation (I5.27p . and o"i/2 times equation (I5.24p . 
This gives 

(5.29) (a 2 j4 - a 2 )g 22 - ^g 2 - b 2 g = -(C(<n) + aC(a 1 ))/2. 

Noting similarity with (I5.20p . we again we consider (I5.29P as a differential equation with 
respect to o 2 in the interval a 2 < erf/ 4. With b ^ 0, t = of/4 — cr 2 , we take 

= a- 2 ) — , 

and we get 

2ty"+(l + (3)y'~2b 2 y = 0. 

Therefore, the solution bounded at the origin is g(a\,a 2 ) = K (ai)X^_iy 2 {b^ a\ / 4 — a 2 ). 
Substituting this back into (15.241) . we get K' = bK, i.e., K(o~i) = e bai , and the resulting 
g{o-i,a 2 ) solves (15.231) . too. 

Thus the general solution bounded at the origin is 



<?K a 2 ) = C x + C 2 e 2b ^ + C 3 e 6CTl Z ( ^ 1)/2 (6 v /(7 1 2 /4 - a 2 ). 

The initial conditions are satisfied when C\ + C 2 + C3 = and 2C 2 + C3 = 1/ (2b 2 ). After 
some calculations, we get the answer. □ 

5.2.2. Gaussian case: a = 0, b = 0. The system of PDEs in Theorem I4.1( ii) simplifies to: 

(5.30) - cr 2 5f 22 = ^g 2 + 1, 

(5.31) 25-12 + o x g 22 = 0. 

Proposition 5.8. Consider the system (15.301) . (15.311) with zero initial conditions. The 
solution is 

(5.32) C<t?/2-(1-C)|<7 2 , 
where C G R is arbitrary. Thus 

(5.33) L(0) = expfc?^-(l-C)|a 2 Y 
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Since the homogeneous system is the same as in the case b 2 — 4a = of Proposition 
I5.6( i). we omit the details. 

We recall the following facts about univariate Laplace transforms, which belong to the 
folklore: 

Proposition B. Suppose L(t) is a Laplace transform of a real random variable X . If L(t) 
is defined on an open interval Q, then L has analytic extension to + iR. If Zq G O + iR 
is a zero of L of order m G IN, then L a cannot be a Laplace transform of a probability 
measure unless ma G IN. 

Proof of Theorem 15.11 Without loss of generality, we may assume that X is non-degenerate, 
and therefore that it is standardized with mean and variance I n . Therefore, ( II. 4p holds 
with some constants a, b, as explained in the introduction, replacing X by —X if needed, 
without loss of generality we may assume b > 0. Therefore the Laplace transform of X 
must be given by one of the six formulas from Proposition I5.6f i-iii). Proposition 15.71 or 
Proposition 15.81 It remains therefore to show that these functions are the Laplace trans- 
forms of probability measures only when the parameters satisfy the conditions listed in 
the corresponding definitions of our Meixner ensembles. The solutions, and our defini- 
tions, depend on the constraints satisfied by parameters a, b, so we need to consider each 
case separately. 

In the proof, we repeatedly compute L{6) on diagonal matrices 

t/2 
±t/2 ' 

(i) Case b 2 > 4a with a < 0. In this case, all potential Laplace transforms for the 
standardized ensembles are given by formula (15.121) . and (15.31) shows that L(6) is 
indeed the Laplace transform if a = —1/(4N) and C\,C^,Cz > 0. It remains to 
show that L{6) is not a Laplace transform if l/(4a) (jL IN or if one of the constants 
is negative. 

Denote a = -l/(4a) > 0. If C 3 = 0, then L(\6j) = (de* + C 2 e~ t ) a , and it 
is a classical fact that a G IN and C\,C% are the binomial probabilities, compare 
Proposition |3]U Therefore, through the remainder of the proof we assume C3 > 0. 
We first show that a G IN; the same reasoning shows that the J0rgensen set of a 
Bernoulli ensemble on is {1, 2, . . . }. We write L(Bj) = (f(s)) a , where 

f(s) = l-C 3 + C 3 Z { ^ 1)/2 (s). 

To show that a G IN, we will apply Proposition [Bl We first observe that (13.241) 
implies 

00 

(5.34) X,(v^) = ^a n ^, 

n=0 

where 

1 

Qn ~ 2 2n n\Y{n + v + 1) ' 
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This shows that T v (^/z) is an entire function of order 

i- lo § n 1 m 
p = hmsup -j = 1/2. 

n-*oc - log a n 

Therefore by a refinement of the little Piccard's theorem [301 Theorem 9.11], the 
equation T u (yfz) = (C 3 — 1)/C 3 has an infinite number of roots. 
We also observe that 

35) -^-l u (y/z) = ^ +1 (v^), 

see (3HI page 479]. 

Coming back to f(s), we see that the function s H- f(s) has an infinite set 
Z of roots. We claim that Z has at least one simple root. If not, Z would be 
included in the set Z x of roots of the derivative f'(s). Now from (15.351) we see 
that Zi is the set of roots of Z(^ +1 )/ 2 which by Lommel's Theorem, lies on the 
imaginary axis ([381 15.25]). From [381 P a g e 199] we see that 

lim T v (it) = 0, 

t— s>±oo 

therefore T v (it) cannot be equal to (C 3 — 1)/C 3 on an infinite subset of Z\. There- 
fore / has at least one simple zero, and from Proposition [B] we deduce that a 
must be an integer. 

The final step is to show that Ci,C 2 ,C 3 > 0. Since the second derivative of 
logX(^r) = « Ml - C 3 + C 3 %_ 1)/2 (t)) is aC 3 /(P + 1), we see that C 3 > 0. 
Next we compute the second derivative of - log L(-9^), which is 



(C 2 C 3 + e 2t C x C z + 4C 1 C 2 e') 



(C 2 + e* (e'd + Cs)) 



2 • 



Since a G IN, the Laplace transform L{-6f) is well defined for all real t. The sign 
of the second derivative for large t is determined by C1C3, showing that C\ > 0. 
The sign of the second derivative as t — > — 00 is determined by C 2 C3, showing 
that C 2 > 0. (Recall that we consider the case C 3 7^ only.) 

Thus Definition 13.21 indeed covers all examples of Meixner ensembles on H 2i( g 
with a < 0. 

Case b 2 > 0, a = 0. In this case, all potential Laplace transforms for the standard- 
ized ensembles are given by formula (I5.26p . From (15 .4p . after passing to standard- 
ized X we confirm that (15. 26j) is indeed a Laplace transform when 1/2 < C < 1, 

To verify that L(8) fails to be a Laplace transform of an ensemble in all other 
cases, we compute 



2 



d_ 



log Lfo/b) 



= 2C-1. 
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Thus we must have C > 1/2. Next, we note that in order for 

^logL (0+/b) = e* (4e*(l - C) + 2C - l) 

to be positive for large t, we must have C < 1. 

This shows that Definition 13.31 covers all examples of Meixner ensembles on 
H 2j/3 with 6 ^ 0, a = 

(iii) Case b 2 > 4a, a > 0. In this case, all potential Laplace transforms for the 
standardized ensembles are given by formula (I5.13p . To see which values of A 
correspond to the Laplace transform (I5.5P of the negative binomial ensemble, we 
take r = l/(4a) and define parameters qi, q 2 and p = 1 — q\ — q 2 as follows 

2 

- = 1- X + b/K, q 2 = (p - l)pA, g! = 1 -p - q t . 
p 

From the admissible range of p given in Remark 15.51 we see that (I5.13P is a Laplace 
transform of an ensemble when 1 < 1 — \ <b/n. 

It remains to show that no other values of A are allowed. The second derivative 

of 

4alogL(0 t -/«) = -log(l- A + AZ (/3 _ 1)/2 (t)) 

at t = is -A/(l + (3), proving that A < 0. 

Next suppose 1 — A > b/n. Then, with p = — A)) < 1 we see that 

(1 — A) cosh t — b/n sinh(t) = (1 — A) (cosh t — p sinh t) 

is an increasing function of t > 0. So L{6f / k) is well defined for all t > 0. We 
now compute the second derivative of 

(5.36) 4a\ogL(6t/K) = - log(A + (1 - A) cosh(t) - 6/Ksinh*). 

We get 

b 2 - k 2 (\ - l) 2 + kX(k{\ - 1) cosh(t) + 6sinh(t)) 
(-kA + k(A - 1) cosh(t) + 6sinh(t)) 2 ' 

Since k(A — 1) + b < 0, this expression fails to be positive for large t, contradicting 
that L is a Laplace transform of a probability measure. Thus we must have 
1 — A < &/«. 

This shows that Definition 13.41 covers all examples of Meixner ensembles on 
H 2)1 a with a > 0, b 2 > 4a. 

(iv) Case a = 0, b = (Proposition 15. 8p . It is clear that Definition 13.51 covers all 
examples of Meixner ensembles on 13-2$ with Laplace transform f)5.33p . 

(v) Case a > 0, b 2 = 4a (Proposition 15. 6K p|). From [28] it follows that Definition 13.61 
covers all examples of Meixner ensembles on H^/? with a > 0, b 2 = 4a. 

(vi) Case b 2 < 4a (Proposition I5.6lflrli|) ). From Proposition 13.141 we see that (I5.14p is 
a Laplace transform of a probability measure when < A < 1, and also when 
A = 1, b = 0. It remains to show that (13.251) fails to be a Laplace transform in 
the remaining cases. 
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Denote a = l/(4a) > 0. The second derivative of 

logL(0 t ") = -«log(l - A + AJ- (/3 _ 1)/2 (0) 

at t = is a\/ (1 + 0), proving that A > 0. 

Next, for A ^ 1 let G (-7r/2,vr/2) be such that sin0 = p/y/(l - A) 2 + p 2 . 
Then the second derivative of 

log L(8+) = -alog(A + (l - A)(cost-p/(l - A) sin*)) 

cos(t + 4 



a log A + (l-A)- 

cos< 



at t = —(f) is 

1-A 

a — 



-1)A + 1 



\ cos t, 

Since A > 0, this shows that A < 1. 

It remains to show that (I5.14p is not a Laplace transform of a probability 
measure when A = 1 and b ^ 0. In this case, consider 

g(s,t) = log£(~(0+ + t -)) = -fe/(4a) -odog (^-i)/ 2 (0 - 6sin(s)/«) . 

At 

Then g(s, t) is well defined on (s, t) G [0, 7r) x (— Ep, £p). We calculate the Hessian 

H(s,t) = det 
at s = and we get 



Using Mathematica, we verify that J~"p_iy 2(^)^(13-1)/ 2(^/3) < 0' so for b 7^ the 
Hessian must change sign over the domain (s, t) G [0, 7r) x ^3). Thus logL(0) 
cannot be a convex function, i.e., when A = 1, b 7^ 0, function L{9) cannot be a 
Laplace transform of a probability measure for any a^0. 

So Definition 13.71 exhausts all the Meixner ensembles on H 2) /3 with b 2 < 4a. 

□ 



6. Additional observations 

6.1. Independence of S and S _1 X 2 S _1 . Bozejko and Bryc [91 Remark 5.8] raise the 
question whether there exists a non-trivial law on positive n x n matrices such that 

(6.1) Z := S _1 X 2 S-\ S := X + Y 

are independent when X and Y are i. i. d. matrices with this law. Such laws could 
provide matrix models for the "free gamma" law. 

Here we answer this question in negative, at least for the laws on positive 2x2 matrices 
which are invariant under orthogonal/unitary/symplectic group. We show that in the case 
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of 2 x 2 positive random matrices, S _1 X 2 S _1 and S are independent only if X arises as 
a gamma random variable multiplied by I2. 

Proposition 6.1. // X, Y G H 2)i g are i.i.d. square-integrable non-degenerate positive 
random matrices such that S = X + Y and S -1 X 2 S -1 are independent and rotation 
invariant, then X = £I 2 and £ has a univariate gamma law. 

Lemma 6.2. Suppose X, Y G H 2)/ g are positive i.i.d. random matrices with the Laplace 
transform 

(6.2) L{9) = (1 +Ati6 + BtT 2 6 + Ctr9 2 )~ p . 

Let S = X + Y. If real random variables det(X)/ det(S) and det S are independent, 
then one of the following cases occurs: 

(i) A = B = C = 0. (Then X = I.J 

(ii) A = ±2\^B, C = 0. (Then X = ^£1 where £ zs univariate gamma with shape 
parameter 2p.) 

(iii) A = ±V2B, C = —B. (Then X = AY where Y has a Wishart law with shape 
parameter p, i.e., with the Laplace transform L{9) = (1 — tr (9) + det8)~ p = 
det (I - 9)- p .) 

Proof. The independence of determinants implies 

E (det X exp tr (OS)) = c E (det S exp tr (6S)) , 

9\ #12 
6*12 #2 



where c = E (f|f ) > 0. With 5 



this becomes 



a 2 a 2 \ w „, / a 2 d 2 



Writing (O as L{6) = l/(f(9)) p , when c ^ ^gri) we § et 

df Of I f df \ 2 ( d 2 f 1 d 2 f 



89x 892 4\89 12 J Cl ^\<90id0 2 4<9fl 2 2 
where 

2cp-l 
Cl 2(2p+l)co-p-l ' 
Solving the resulting equations for A, B, C, c\ we verify that the only non-zero solutions 
correspond either to c x = 2, A = ±2^B, C = with f(9) = (1 + Atr (9)/2) 2 , or to 
d = 2/3 with A = ±y/2B, C = -B. 

When c = 2( g^ 1} , the equations give C = 2B so /(0) = 1 + Atr6 + B(ti 2 9 + 

2tr# 2 ). Using [281 Theorem 3.1] one can verify that in this case L{9) is not a Laplace 
transform. □ 

Proof of Proposition l6A[ Denote Z = S -1 X 2 S -1 . Note that this is a positive-define ma- 
trix. First we note that for positive matrices the Laplace transform is defined on 9 G H 
when —9 is positive. 
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Next we note that rotation invariance of the law of X implies that E(X) = al with 
a > 0, E(X 2 ) = (a 2 + (3 2 )I with (3 2 > 0, and E(Z) = c I. Therefore, 

E(X 2 |S) = E(SZS|S) = c S 2 . 

Taking the expected value of both sides, cq = 2 {p i +2a 2 ) ■ ^ e a ^ so no ^ e ^ na ^ by exchange- 
ability, E(X 2 |S) = E(Y 2 |S). So 

E((X - Y) 2 |S) = E(2X 2 + 2Y 2 - S 2 |S) = (4c - 1)S 2 . 

Passing to standardized matrices X' = (X — aI)//3, we see that (11.4j) holds with o 2 = 
4a = a 2 /(3 2 . Thus by Proposition I5.6( i). 

E (exp tr (0X)) = e atr9 E (exp tr (/30X')) 

and equation (I6.2p holds with p = (3 2 /a 2 > 1. 

Since det(X)/ det(S) = VaetZ is independent of det S, by Lemma |6T2| there are only 
two non-degenerate choices for the law of X. It is known that (ll.4|) does not hold for the 
Wishart law, see [271 PS- 582], which excludes the law from Lemma [6.2( iii). Thus X is 
the identity matrix multiplied by a Gamma random variable. □ 

6.2. Some series solutions of the general system of PDEs. The Laplace transform 
of Bernoulli ensembles can be computed readily. For example, if P is the orthogonal 
projection onto a randomly rotated line through the origin in C 3 then from [26J we have 
for 9 = Diag[9i,e 2 ,e 3 ), 

(6.3) Eexp(0|P) = 2— — — — + 2- 



^-.72,^2-03) (01-03) (03 -02) 

However, it does not seem to be obvious how to express this Laplace transform in terms of 
the elementary symmetric functions o\ = 01 + 02 + 03, cr 2 = 6i6 2 + 0i03 + 2 03, 03 = 010203- 
The following result gives the expansion for rank 1 projections in H n / g for all n > 2. When 
combined with properties of Meixner ensembles from Section[3J this can be used to produce 
some solutions of the system of PDEs, see e.g., Corollary 16.41 
Recall the Pochhammer symbol 

(6.4) (b) k = ^±^ = b(b+l)...(b + k-l). 

Proposition 6.3. If Pi is an random projection o/H n)/ g invariant by rotation with trace 

1, then Ee^ F ^ = L n (9) where 

(6.5) 

OO / i\fc ( — 1)" l+ V 2+f3+- fM") 

L ^ = E Wf E — „.„ ... V1+VM - °ne)°?{e)°? (0) • • • 

k=0 V 2 Ik yi+2y 2 +3i/ 3 +---=fc 1- 2 ' 3 ' ' ' ' 

Since P2 = I3 — Pi for n = 3, we have the following. 

Corollary 6.4. The Laplace transform of the Bernoulli ensemble on H 3j/ 3 with parameters 
91,92,93 ^ 

L{9) = q + qiL 3 (9) + +q 2 e^L 3 (-9) + q 3 e^ . 
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For the proof of the proposition we use facts about Jack polynomials (spherical poly- 
nomials) which require additional notation. 

6.2.1. Partitions. A partition k = (ki, k 2 , ■ ■ ■ , K) of k into at most r parts is the sequence 
k\ > k 2 > ■ ■ • > k r > such that k\ + k 2 + ■ — h k r = k. We write k h k. Other standard 
notation is \k\ = k± + k 2 + • • ■ (which is just k). We denote the length (the number of 
parts) of a partition k by £{k) = #{j : kj > 0}. 

6.2.2. Pochhammer symbols. For a partition A h fc, we need the /3-Pochhammer symbol: 

(6.6) (^^rii 6 -^- 1 ' 



6.2.3. JacA; polynomials. For a not necessarily symmetric n x n matrix 9 with eigen- 
values 9i,9 2 , . . . ,9 n , the Jack polynomials {J K (0;/3)} are the symmetric polynomials in 
9i, 9 2 , . . . , 9 n , indexed by the partitions k of k = \k\ into r = £(k) < n parts. In [2T?1 
formula (1)] the Jack polynomial J K (#; /3) is defined as the coefficient at t±t 2 . . .t r in the 
expansion of 

(2//3) r det(J n - h9 kl - t 2 9 k2 t r 9 kr )~^ 2 . 

It is convenient to set J K (0; 0) — if the partition k has more than n parts. (The usual 
notation for J K (x; (3) is C K (x; a) with a = 2/(3.) 
We need several properties of Jack polynomials. 

Proposition C. (i) If £(k) < n and A,B are n x n matrices, then 
(6.7) / UUAU>B;WU- UA - mP) 



where the integral is over the normalized Haar measure on the orthogonal group 
O n ((3 = I), on the unitary group U n ((3 = 2), or on the symplectic group K n 

(ii) If A is an n x n matrix and k > 0, then 
(6.8) (Ktr(A)) fc = ^J K (A;/3). 



(6.9) 



iii) If £(X) < min{m, n}, then 

$\{lm,n',(3) \~2~)\;t 



Ja(I»;0 (f) A;/3 

(If£{\) > m then J A (I m , n ; (3) = O.J 
iv) For a one-part partition A = (k) h k, we have 



V2/fc i/i+2i/ 2 +3^3+---=fc ' - 
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Formula (16. 7ft appears in [18], where it is attributed to [HI (45)], see also p2J Corollary 
XI. 3. 2]. Formula f !6.8p appears in [33 Proposition 2.3], see also [HH page 234]. Ref. [221 
(21)] gives J K (I m ; 1). Formula (KTTOjl is [33 Proposition 2.2(c)]. 

We are now ready to prove Proposition 16.31 

Proof of Proposition 1 6. 3[ Recall that Pi = U*li >n U, where Ii )7l is the identity matrix 
padded with zeros to make it into the n x n matrix. From (16.81) . 



oo 1 00 1 

Eexp^jPi) = ^-Etr fc (£W*I lin ) = ^ - ^ $ K (WU*I hn ; (3) 
By (16. 7p . this becomes 



k=0 k=0 fthfc 



7J 



J K (0;0)JU(Ii, n ;0) 



Since J K (Ii, n ; /3) = when > 1, this gives 

1 J( fc )(^;/3)J( fc) (Ii, n ;/3) 



oo 



Eexp^P^V 

^fc! J(fc)(In;/3) 



Now by (16.91) this becomes 

oo rp\ 

Eexp(^|P 1 ) = ^-^-J (fc) (^/3), 

fe=Q \ 2 Sk 

and formula ( 16. 10}) gives the answer. 

□ 

6.3. Connection with free probability. The free probability analog of Meixner family 
was introduced by Anshelevich, Saitoh and Yoshida [31 EE] as the orthogonality measure of 
certain orthogonal polynomials. These free Meixner laws include the Kesten law (a special 
case of the free binomial law, which is the free additive convolution of the Bernoulli law), 
the Marchenko-Pastur law (also called the free-Poisson law), and the Wigner semi-circle 
law (a free probability analog of the normal law). The free Meixner family of laws again 
appeared as the laws characterized by a quadratic regression property in free probability 
[9], in a class of Markov processes, and as generating measures of Cauchy-Stieltjes kernel 
families with quadratic variance function [10]. See also [5], [6] and [HI Theorem 4.3]. 

Except for the free binomial law, the remaining free-Meixner laws are infinitely divisible 
with respect to free additive convolution and appear as limit laws of large dimensional 
random matrices [3 [12]- However, this correspondence is based on the Bercovici-Pata 
bijection [8], which, as observed by Anshelevich [31 page 241] is different from the corre- 
spondence based on the kernel families or regression properties. 

6.4. Connection with Euclidean Jordan algebras. The natural framework of the 
Anshelevich question is rather the one of the Euclidean Jordan algebras well described 
in the book of Faraut and Koranyi [16]. The spaces H nj( g for /3 = 1,2,4 are the 3 more 
important types of (irreducible) Euclidean Jordan algebras. But there are two more types 
to be considered: the exceptional 27 dimensional Albert algebra which can be roughly be 
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X\ V 
t V X 2 



considered as a space of (3,3) Hermitian matrices on octonions, and the Lorentz algebra 
on R 2 x E where E is a Euclidean space of dimension d — 2. In this algebra, the product 
of (x 1 ,x 2 ,v) and (y 1 ,y 2 ,w) is 

(xxyi + (v,w),x 2 y 2 + (v,w),-(yi+y 2 )v+ ~(xi + x 2 )w) 

which leads to the square of (x%, x 2 , v) as (x\ + \\v\\ 2 , x\ + \v || 2 , 2{x\ + x 2 )v). A good way 

to memorize this product is to write formally (xi, x 2 , v) as a 2 x 2 matrix X - 

doing the same with Y for (yi, y 2 , w) and to consider the Jordan product 

X-Y = -(XY + YX). 

If d = (3 + 2 with (3 = 1,2,4 this algebra is H 2/ 3. For a Jordan algebra V, rotational 
invariance in H 2 ./3 is generalized to the invariance by a certain compact group K acting 
on V, as described in Faraut Koranyi [T6"l page 6 and 55]. The Anshelevich problem can 
be solved for the Lorentz algebra for any d > 1 since the rank is 2, and the preceding 
study for H2./3 extends easily to this case. Since the consideration of the Lorentz and 
Albert algebras can be seen as rather academic, we have refrained to place our paper in 
this framework. 

6.5. Unresolved questions. 

6.1. For fixed A, B, C the research of the set M(A, B, C) of the ensembles on M. n ^ which 
satisfies E((X — Y) 2 |S) = AS 2 + BS + CI n led for n — 1 to the Laha Lukacs study 
which shows M(A, B, C) is a natural exponential family. For n > 2 the structure of the 
M(A, B, C) is not really understood: is it possible to pass easily from the knowledge 
of one element of M(A, B, C) to a knowledge of the whole set as it is the case for an 
exponential family? The answer is hidden in the structure of the PDE system. 

6.2. The dimension of the set of solutions of the PDE systems from Theorem 14. II is equal 
4 if n = 2 and is unknown if n > 3. The subset of solutions that are analytic at zero has 
dimension 3 for n = 2, and is unknown for n > 3. A related system of PDEs in [21], has 
a one-dimensional set of solutions analytic at 0. 

6.3. Since the free binomial law is well defined for all real values of iV > 1, it would be 
interesting to determine the J0rgensen set of the Bernoulli ensembles for n > 2. 

6.4. For n > 3, it is not known whether a Meixner ensemble with parameters A = 
-1/(2N - 1), B = 2N/(2N - 1), C = is a binomial ensemble from Definition ET21 
Similarly, the converses to Propositions 13.61 and 13.71 are not known beyond dimension 
n = 2 which is known from Theorem 15.11 

6.5. For n > 3, do we have non-Gaussian ensembles on such that if X, Y are i.i.d. 
then E((X - Y) 2 |X + Y) = 2I n ? 

6.6. The distribution of eigenvalues as the dimension n — > 00 was studied in [12] for a 
special case of Poisson ensembles. It is interesting to study this limit for other Meixner 
ensembles. 
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